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SERIES REPRESENTATION OF POWER FUNCTION 


KOLOSOV PETRO 


Abstract. In this paper described numerical expansion of natural-valued 
power function x n , in point x = xo, where n,x o - positive integers. Applying 
numerical methods, the calculus of finite differences, particular pattern, that 
is sequence A287326 in OEIS, which shows the expansion of perfect cube n as 
row sum over /c, 0 < k < n — lis generalized, obtained results are applied to 
show expansion of monomial 7l 2m + 1 , m = 0,1,2, ...,N. Additionally, relation 
between Faulhaber’s sum J2 77,771 an d finite differences of power are shown in 
section 4. 
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1. Introduction and main results 

In this paper particular pattern, that is triangle A287326 in OEIS, m , which 
shows necessary items to expand perfect cube n as row sum is generalized and 
obtained results are applied on expansion of monomial /(n) = n m , ( n,m ) £ N. 
The coefficient M\{n,k) is fc-th item of n-th row of triangle A287326, such that 
sum of the n-th row of A287326 over 0 < k < n— lis perfect cube. First, let review 
and basically describe Newton’s Binomial Theorem, since our coefficient M\(n,k) 
is derived from finite difference of perfect cubes, which is taken regarding Binomial 
expansion. In elementary algebra, the Binomial theorem describes the algebraic 
expansion of powers of a binomial. The theorem describes expanding of the power 
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of (x + y) n into a sum involving terms of the form ax b y c where the exponents b 
and c are nonnegative integers with b + c = n, and the coefficient a of each term is 
a specific positive integer depending on n and b. The coefficient a in the term of 
ax b y c is known as the Binomial coefficient. The main properties of the Binomial 
Theorem are next 


Properties 1.1. Binomial Theorem properties 

(1) The powers of x go down until it reaches Xq = 1 starting value is n (the n 
in (x + y) n ) 

(2) The powers of y go up from 0 (y° = 1) until it reaches n (also n in (x + y) n ) 

(3) The n-th row of the Pascal’s Triangle (see [I], |12j ) will be the coefficients 
of the expanded binomial. 

(4) For each line, the number of products (i.e. the sum of the coefficients) is 
equal to x + 1 

(5) For each line, the number of product groups is equal to 2 n 

According to the Binomial theorem, it is possible to expand any power of x + y 
into a sum involving Binomial coefficients 

(1.2) (x + y) n = Z! (?)*"" V 

fc=o ' ' 

Let expand monomial /( x) = x n , where x and n are positive integers, applying a 
finite difference operator 

Lemma 1.3. Power function could be represented as discrete integral of its first 
order finite difference 

jkh n ~ 1 +h n 

-•---' 

Ah[x n ]=(x+h) n —x n 

x —1 n / \ 

= EE x,neN, heR 

j =0 k=1 ' ' 

Or, by means of Fundamental Theorem of Calculus 

X-l pk+l 35 — 1 

nt n ~ 1 dt = J2 / nt n ~ l dt = Y^(k + l) n -k n 

fc=o fc=o 

Let describe the derivation of the sequence A287326 in OEIS, which shows the 
expansion of perfect cube n as row sum over k , 0 < k < n — 1. First, review a 



i —n \ / \ 
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difference table of perfect cubes (d, eq. 7) 


(1.4) 


n 

A u (n 3 ) 

A i (n 3 ) 

A 2 (n a ) 

A 3 (n 3 ) 

0 

0 

1 

6 

6 

1 

1 

7 

12 

6 

2 

8 

19 

18 

6 

3 

27 

37 

24 

6 

4 

64 

61 

30 

6 

5 

125 

91 

36 

6 

6 

216 

127 

42 

6 

7 

343 

169 

48 

6 

8 

512 

217 

54 


9 

729 

271 



10 

1000 





Table 1: Difference table of perfect cubes n, 0 < n < 10 of order k, 0 < k < 3. 
Note that increment h is set to be h = 1 and k > 2-order difference is taken 
regarding to 0 , Hi- Reviewing |F igure (1.4) | we have noticed that 

(1.5) A(0 3 ) = 1 + 31-0 

A(l 3 ) = 1 + 3!-0 + 31-1 

A(2 3 ) = 1 + 3!-0 + 3!-1 + 3!-2 

A(3 3 ) = 1 + 3!-0 + 3!-1 + 3!-2 + 3!-3 


A(n 3 ) = 1 + 3! • 0 + 3! • 1 + 3! • 2 + 3! • 3 H-+ 3! • n 

Also, another interesting thing was observed, according to Donald Knuth, the per¬ 
fect cube n is 

f n + 1 N 


( 1 . 6 ) 


n — 6 


+ 


3 ) \1 

Returning to expression m, we can find that 

Ai + 1 


(1.7) 
Then 

( 1 . 8 ) 


A(n 3 ) = 6 


{A 2 (?r 3 )}„ =0 = { 6 



= {6,12,18,24,..., 6f} 


—0 ' n —0 


which fit according to A 2 (n 3 ) values from Figure (1.4) Therefore, we represent 
perfect cube n as 

(1.9) n 3 = (1 + 3! ■ 0) + (1 + 3! • 0 + 3! • 1) + • • • + (1 + 3! ■ 0 + • • • + 3! • (n - 1)) 
Generalizing above expression, we have 

(1.10) n 3 = n + {n — 0) • 3! • 0 + (n — 1) • 3! • 1 + • • • + {n — {n — 1)) • 3! • (n — 1) 
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Provided that n is natural. Now we apply a compact sigma notation on (11.101) . 
hereby 

n— 1 

(1.11) n 3 = n + 3! • nk — 3! • k 2 

k -o 

As sum 3! • nk — 3! • k 2 consists of n terms, we have right to move n from 

(11.111) under sigma notation, we get 

n— 1 

(1.12) n 3 = ^ 3! • nfc — 3! • fc 2 + 1 

fc=o 


Property 1.13. Let be a sets A(n ) := {1, 2, ..., n}, B{n ) := {0, 1, ..., n}, 
C(n ) := (0, 1, ..., n — 1}, let be expression \1.1‘2\) defined as 

T{n , C(n)) := ^ 3! • nk - 3! ■ fc 2 + 1 

feGC(n) 

where x is natural-valued variable and C{n) is iteration set of I11.12\) . then we have 
equality 

(1.14) T(n, A{n)) = T{n, C(n)) 

Let review and define expression \1.10\) as 


then 


U(n, C(n)) := n + 3! • ^ nk — k 2 

fcGC(n) 


(1.15) U(n, A(n))=U(n, B{n)) = U{n, C{n)) 

Other words, changing of iteration sets of i 1.1(A) and \1.12\) by A(n ), B(n), C(n) 
and A(n), C(n), respectively, doesn’t change resulting value for each natural x. 

Proof. Let be a plot y(n , k) = 3! ■ nk — 3! • k 2 + 1, tgR, 0 < k < 10, given n = 10 



Figure 2. Plot of y(n, k) = 3! • nk — 3! • k 2 + 1, k G R, 0 < k < n , where n = 10. 
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Obviously, being a parabolic function, it’s symmetrical over hence equivalent 
T(n, A(n)) = T(n, C{n)) follows. Reviewing (11.101) and denote u(n,t) = n 2 t — 
t 2 , we can conclude, that u(n, 0) = u(n,n ) = 0, then equality of U(n, A(n)) = 
U(n, B(n )) = U(n, C(n)) immediately follows. This completes the proof. □ 

Review above property (11.131) . Let be an example of triangle built using 


(1.16) 


y(n , k) = 3! • nk — 3! • k 2 + 1, 0 < k < 


over n from 0 to n = 4, where n denotes corresponding row and k shows the item 
of row n. 


Row 0 
Row 1 
Row 2 
Row 3 
Row 4 


1 


1 

1 1 

17 1 

13 13 

19 25 19 


1 


(1.17) 

Figure 3. Triangle generated by (11.161) from 0 to n = 4, sequence A287326 in 

OEIS, HU- 

Note that n-th row sum of Triangle (11.171) over 0 < k < n — 1 returns perfect 
cube n. We can see that each row with respect to variable n = 0, 1, 2, 3, 4,..., 
has Binomial distribution of row terms. One could compare Triangle (11.171) with 
Pascal’s triangle in. m 


Row 0 
Row 1 
Row 2 
Row 3 
Row 4 


1 


1 


1 4 6 4 1 

Figure 4. Pascal’s triangle up to forth row, sequence A007318 in OEIS. 

1.1. Properties of Triangle (1.17) and other expansions. Review the triangle 
(flTTl) . define the fc-th, 0 < k < n, item of n-th row of triangle as 

Definition 1.18. 


(1.19) 


Mi(n, k ) := 3! • nk — 3! • k 2 + 1 = 6(n — k)k + 1, 0 < k < 


Let us approach to show a few properties of triangle (11.171) and Mi(n, k). 

Properties 1.20. Properties of triangle |_?.17| h 

(1) Summation of items M\{n, k) of n-th row of triangle {1.1 7| ) over k from 0 
to n — 1 returns perfect cube n as binomial of the form 

n— 1 

y, Afi(n, k) = A 0 , n n - B 0 ^ n = n 3 , n > 0 

k =0 


( 1 . 21 ) 
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Since the property |T7 fj) holds, \1.Z1\) could be rewritten as 

n 

(1.22) M\ ( n , k) = Ai tU n — Bi }U = n 3 , n > 0 

fc=l 

where A^-j n and B^j n - integers depending on variable n £ N and on 
sets U(n), S(n), respectively. Note that Ay^n ^ Ai tH n, Bq , n ^ B i, n 

(2) Relation between Ao, n and A\^ n 

Ao,n+i = A\ n , n > 1 

(3) Summation of items Mi(n,k) of n-th row of triangle \1.1 7| ) over k from 0 
to n returns n 3 + 1 

n 

(1.23) ^Mi(n,ft)=n 3 + 1 

fc—o 

(4) Mrst itom of each row’s number corresponding to central polygonal numbers 
sequence a(n) = " + 2 n+2 (sequence A000124 in OEIS, [T3]^ returns finite 
difference of consequent perfect cubes. For example, let be a k-th row of 
triangle \1.11\ , such that k = - + 2 ra+2 , n = 0,1, 2,..., t/ien item 

(1.24) Mi tZLlti^ = (n + l) 3 — n 3 , n> 0 

(5) Items of |1.17| ) have Binomial distribution over rows. 

(6) The linear recurrence, for any k and n > 0 

(1.25) 2Mi(n, fc) = Mi(n +1, fc) + Mi(n — 1, fc) 

(ZTiis linear recurrence is direct result of second order binomial transform of 
M\(n,k) over n. 

(7) Linear recurrence, for each n > k 

(1.26) 2Mi(n, k ) = M x (2n - fc, fc) + Mi(2n - fc, 0) 

(8) From (TZfj ) follows that 

n —1 n —1 / 2 I I o \ 

(1.27) ^ Mi(n,k) = ^2 M i (--7T-,1 j = ™ 3 , '« > 0 

fc=0 fc=o ' ' 

(9) Triangle dl,17| ) *s symmetric, i.e 

(1.28) Mi(n, fc) = M\(n, n — fc) 

As its noticed in (11.211) . summation of each n-th row of Triangle (11.171) from 0 
to n — 1 returns perfect cube n, then, by properties (11.241) . (11.251) . (11.261) . for each 
positive integer n the n m , m = 0,1, 2,... could be found via multiplication of each 
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term of (11.2111 by ri 


m —3 


(1.29 yi m = y M\{n, k)n m 3 = - y] [M\(n + 1, k) + M\{n — 1, k )] ri 


m —3 


k=0 
n —1 


fc=0 


A 1 

E - [Mi(2n - k, k ) + Afi(2n - k, 0)] n m ~ 3 


k—0 
n —1 


= E> 


1 _ _ / n 2 + n + 2 


/c—0 

n— 1 _ r- / o 

n + n 


■ 1 ) n 


m —3 


E- 

^ 2 

k—0 

n —1 

E- 
^ 2 

k-0 


Mi 


Mi 


, 1 + Mi 


n 2 + ri + 4 


_ m—3 


m + 1 
2 


, 1 + Mr 


m + 1 
2 


,1 


„m— 3 


To show other representation of monomial ro m , m = 0,1,2,..., review (11.111) . let 
move n in n + 3! y\ — A: 2 under the sum operator and change iteration set from 
{0, n — 1} to {1, n — 1}, then we get 

n— 1 


(1.30) 


n 


= ^ 3! • nk — 3! • k 2 


fc =i 


(n - 1) 


, «/l 


Review right part of (11.301) . let the term be written as 1 = l + 

given the power m > 3, multiplying each term of (11.301) by n m ~ 3 we can observe 
that 

n —1 

(1.31) n m - 1 = E Ml ( n ’ + « m_4 + n m ~ 5 + • • • + n + 1 

k= 1 

Applying properties (ll.24[) , (11.25p , (|1.26p , we can rewrite (|1.3ip as 

n— 1 


n m - 1 = 


— [Mi (2n — k, k) + Mi (2n — /c, 0)] n m 3 4- n m 4 + • • • + n + 1 
k =l z 

n— 1 y 

(1.32) = i [Mi(n + 1, fc) + Mi(n — 1, k)] n m ~ 3 + n m ~ 4 + ■■•-)- n +1 

k =1 Z 
n—1 1 

E 2 [ Ml 0 + Ml ("EfE x ) 


n"* -3 + n m ~ 4 + • • • + n + 1 


m—0 
n—1 


E m >( 


^+”+ 2 ,1) n m ~ 3 + n m ~ 4 + ••• + ! 


fc=i 


Review (11.311) . let move 1 from left part of (11.311) under sum operator in the right 
part, therefore we add a term to initial function Mi(n, k)n m ~ 3 +n m ~ 4 +n m ~ 5 + 

- in + 1. By means of expansion — = 1 +n + n 2 +n 3 H-, the (11.311) 

could be rewritten respectively 

n— 1 

n m = J2 M ^n, m)n m ~ 3 + n m ~ 4 + -- - + n+ l- l-n 2 -n 3 - 

m= 1 


(1.33) 
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Particularizing (11,3311 we have 

n— 1 

(1.34) n m = ^ [Mi(n, m) — 1] 

' n 2 +n + 2 


3 2 _ra—1 

n — n — 71 


m—1 
n—1 


= E 


Mi 


ra—1 L 

Hereby, the follow question is stated 


1-1 


n m — 3 - n m_2 - n” 1 " 1 


Question 1.35. Has the Triangle anj/ analogs in order to receive monomial 

x* t > 3 as row sum? Is it exist M v (n, k), v ^ 1, such that 

n —1 

M v (n, k) = n*, v j-tl 

k=0 


2. Generalization of sequence A287326 


Considering the OEIS sequences A300656 and A300785, dquestion 1.31| ) can be 
answered positively, since the sequences A300656 and A300785 are analogs of Tri¬ 
angle (11.171) . that show fifth and seventh powers as row sums. To generalize our 
sequence A287326 for each odd power 2 m +1, m = 0,1, 2 ... we have to review the 
generating formulas of sequences A287326, A300656, A300785 as 

n— 1 m 

(2.1) A m j(n — k) j k j = n 2m+1 , m = 1,2,3 

k—0 j =0 

Where A m j is unknown. For example, generating functions of our sequences 
A287326, A300656, A300785 are 

f 6k(n -k) + 1, for A287326 

(2.2) < 30fc 2 (n — k) 2 + 1, for A300656 

[ 140fc 3 (n - k) 3 - 14fc(n - k) + 1, for A300785 

Reviewing (12.21) we observe that coefficients are {6,1}, {30,0,1}, {140,0, —14,1} in 
each generating function of A287326, A300656, A300785, respectively. To generalize 
above results over odd powers 2m +1, m = 0,1, 2 ... one has to solve the system of 
equations. Define the part of (O) as 


(2.3) M m (n,k ) 

Hereby let be a system 


A m j(n — k) j k j , m > 0 

j—0 


(2.4) 


'M m (l, 0) = l 2m+1 
M m (2,0) + M m (2,l) = 2 2m+1 
< M m ( 3,0) + M m ( 3,1) + M m ( 3,2) = 3 2m+1 


M m (t, 0) + M m (t, 1) H-1- M m (t,t- 1) = t 2m+1 , t > m 
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Solving system m we receive a set of coefficients_{A mi Oj • ••, An,m} for each 
particular t, such that m holds. List of solutions^] of system (Efl) is split and 
assigned to OEIS under the numbers A302971 (numerators of A m j) and A304042 
(denominators of A m j). To reach recurrent formula of A m j, first let fix the unused 
values A m j = 0, for j < 0 or j > m, so we don’t need to care about the summation 
range for j, then by expanding (n — k) J and using Faulhaber’s formula, we get 

n— 1 n —1 oo 

(2.5) £)(n- kyk 3 = 

k =0 i 
i-i (- 1 ) i 


k =0 


i 

oo 

= £ 


n i_i (-l) < fe i+i 


i,t 


i+j + l 

7^TtC + ) + 1 )^ +1 --£(( 


Ml 

« + j + 1" 


■/j 


*+j+l 


1.J-* 


where St are Bernoulli numbers [14] . Now, we notice that 


( 2 . 6 ) 


£ 


(-1)' 


i J i + j + 1 


i+j + l 
t 


= J A+^1’ iH -° ; 

— ( 27 -- t+1 )> if<>0 


t V2j-t+lj 

In particular, the last sum is zero for 0 < t < j. To combine or cancel identical 
terms across the two sums from right part of (12.511 more easily we introducing 
£ = 2 j + 1 — t and £ = j — i, respectively, we get 

n— 1 


(2.7) 


— kyk j = 

k =0 

oo 

- £ 


i 


£=—c 


(2j + l)Cj) 

(-1 m 


,2j+! 


£ 

£=—oo 


Ml 


2j + i-e 


70 i 


B2j+i~en 


1 


(2j + l)( 2 /) 


2j + l-t 

OO 

«+! + 2 £ 


n 


Ml 


-J I odd l 

Note that binomial coefficient is defined as 


2j +1 — t \e 


e 


k\(n—k)\ > if 0 < k < n; 


I 0, otherwise 

Therefore, we don’t use strict limits on sums of above derivation as it’s much easier 
to review each sum as summing from —oo to +oo (unless specified otherwise), where 
only a finite number of terms are nonzero. Now, using the definition of A m j we 
obtain the following identity for polynomials in n 


(2-8) £1 


1 


(2j + l)(£ 


,2j+i 


2 £] A rrl ] 

j , odd Z 


(-i i 


£ J 2j + 1 - £ 


B2j+i-in 


= n 


2m+l 


^One can produce a list of solutions of system (2.4) up to t = 11 using Mathematica code 
solutions_system_2_4.txt, 1241 . 
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Taking the coefficient of n 2m+1 in above expression, we get A m ^ m = (2 m + 1)( 2 ™), 
and taking the coefficient of x 2d+1 for an integer d in the range m/2 < d < m we 
get An.d = 0. Taking the coefficient of n 2d+1 in (2.8) for m/4 < d < m/2 , we get 


(2.9) 


A 


1 

(2d+l)( 2 d d ) 


+ 2(2 m + 1) 



/ m \ (-l) m 
\2d +1/ 2m — 2d 


i.e 

( 2 . 10 ) 


An,<2 = (-l) r 


(2 m + 1)! 


d\d\m\(m — 2d — 1)! m — d 


II 


2m—2d- 


Continue similarly, we can express An,, for each integer j in range m/2 s+1 < 
j < m/2 s (iterating consecutively s = 1,2,...) via previously determined values of 
An,< 2 , d < j as follows 


(2.ii) a„,, = (2 j+ i)Q) ^ ^ ^( 2 j d + i) 

The same formula holds also for d = 0. Hereby, we have shown that EH) holds for 
each m > 0. 


Definition 2.12. We define here a generalized sequence of coefficients An,,, such 
that Efe=o J2jLo A rn,j( n ~ kyk J = n 2m+1 , n > 0, m = 0,1,2,... 

{ 0, if j < 0 or j > m 

(A + 1)( j) Dd=2,+1 Am,d{ 2 j+i)^ d-j -®2d-2,, if 0 < j <m 
( 2 J + !)( 2 /), if j = m 


First five rows of triangle generated by A m j are 


1 


(2.13) 


1 6 

1 0 30 

1 -14 0 140 

1 -120 0 0 630 

1 -1386 660 0 0 2772 


Figure 5. Triangle generated by A m j, 0 < j < m. 

Note that starting from row m > 11 the terms of Triangle (12.131) consist rational 
numbers, for example, Hup = 800361655623.60, therefore it’s split into two se¬ 
quences A302971, A304042 that show numerators and denominators of A m j, 0 < 
j < m , respectively. To verify the terms that definition (12.121) produces one should 
refer to Mathematica cod(0. 


Mef_2_12.txt [25] 
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2.1. Properties of M m (n,k) and A m j. Here we show a few properties of defini¬ 


tion M m (n , k), some of them correlates with properties of Mi (n, k) in subsection (1.1) 

(1) Sum of A m j, m > 0 gives 


(2.14) 


= 2 2m+1 — 1 


(2) Generalization of (12.11) for all positive integers m > 0 


n— 1 


~ k Y k3 = 


2ra+l 


k=0 j 

Can be verified via Mathematica codt0. 

(3) Similarly to particular property (11.281) . items of {M m (ri, fc)}£ =0 , m > 0 is 
symmetric, i.e 

(2.15) M m (n,k) = M m {n,n- k) 

(4) From (12.151) immediately follows 

n—1 n 

(2.16) M m (n, k ) = ^ M m (n, k ) = n 2m+1 , n > 0, m > 0 

fc =0 fc =1 

(5) Generalization of linear recurrence (11.251) . that is 2M\ (n. k) = Mi(n + 
1, fc) + M\(n — 1, fc) 

t / t -X- I \ 

M m (n + t — k, k) = 0, n > 0, t > m 


(2.17) £(-l) 

k—0 

then 
+ 1 
t 


t + 1 
k 


(2.18) 


M m (n,k) = 


EU(-l) fe+1 (T)^n> + t - fc, fc), f = even 

El=o(- 1 ) fe (T) M ™( n + f ~ *,&)> t = odd 

(6) j4 m>m , m = 0,1, 2,... are terms of A002457. 

2.2. Generalized Binomial Series by means of properties (1.21), (1.22). 

Reviewing properties (11.211) and (11.221) . we can say that for each natural n holds 

( 2 - 19 ) = A^n™- 2 - B^ n n m ~ 3 

Note that A^j n and A m j are different definitions. Rewrite the right part of (12.191) 
regarding to itself as recursion 

n m = - ^dTT^^ -5 ) - 

= A hn nm ~ 4 ~ 2 Aop , n % I ,^ 5 + K I ,,.""*" 6 

Reviewing above expression we can observe Binomial coefficients before each A^j n ■ 
Bqj n . Continuous j-times recursion of right part of (12.191) gives us 


( 2 . 20 ) 


= S' - 1 > t (i) A ia h n B h,^ m ^ t ’ i a 0 


k—0 


Solutions Aqj n , Bqj n of equation (|2.19[) are listed in follow table 


c expression_2_l. txt, m■ 
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X 

Ao f x 



Bi,x 

1 

i 

0 

6 

5 

2 

6 

4 

18 

28 

3 

18 

25 

36 

81 

4 

36 

80 

60 

176 

5 

60 

175 

90 

325 

6 

90 

324 

126 

540 

7 

126 

539 

168 

833 

8 

168 

832 

216 

1216 

9 

216 

1215 

270 

1701 

10 

270 

1700 

330 

2300 


Table 9. Array of coefficients A^j n , B^j n given n = 1,10. 

Sequence Ai >x is generated by 3n 2 + 3n, sequence A028896 in OEIS, |23]. Sequence 
_Bi, n is generated by 2n 3 + 3n 2 , sequence A275709 in OEIS, [2D], 


3. Relation between Pascal’s Triangle and Hypercubes 


In this section let review and generalize well known fact about connection be¬ 
tween row sums of Pascal triangle and 2—dimension Hypercube, recall property 


(3.1) 


E 

fc=o 


= 2 r ‘ 


Now, let multiply each k-th term of of n-th row of Pascal’s triangle [I] by 2 fc 

1 

1 2 
1 4 4 

1 6 12 8 
1 8 24 32 16 

Figure 10. Triangle built by (^) • 2 fe , 0 < k < n < 4. 

We can notice that 


(3-2) E 

k=0 

Hereby, let be theorem 


2 k = 3 n , 0 < k < n, (n, k) e N 


Theorem 3.3. Volume of n-dimension hypercube with length m could be calculated 
as 

n k 


(3-4) 


m 


= ££ wu 

k=0j=0 v 7 


(- 1 ) 


k ~ j m j 


where m and n - positive integers. 

Proof. Recall induction over m, in m is shown a well-known example for m = 2. 


(3.5) 


k—0 


2 -=e : ( 2 - 1 ) 
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Review (13.51) and suppose that 

(3.6) (2_+l)" = Ef”V.( 2 - 1 ) + 1 ) J 

k—0 x ' S v / 


m—3 


m— 1 


And, obviously, this statement holds by means of Newton’s Binomial Theorem [2j, 
[3] given to = 3, more detailed, recall expansion for (x + l) n to show it. 


(3.7) 


(* + !)" = £ 


k—0 


Substituting x = 2 to (13771) we have reached m- 

Next, let show example for each to G N. Recall Binomial theorem to show this 


(3.8) to" = £P)(to- l) k 

k =0 x ' 

Hereby, for m + 1 we receive Binomial theorem again 

n 

(3.9) (m + 1)" = 


k—0 


Review result from (13.81) and substituting Binomial expansion Yj=o (j)(~ 0™ km7 
instead (to — l) fc we receive desired result 


(3.10) m n = Y, 


( to . — 1 )* 


k =0 x / , /u , k—0 


k ~irr,i 


=e : e 




(- 1 ) 




= EE (-D 

A:=0j=0 v 7 v// 

This completes the proof. □ 

The (13.51) is analog of MacMillan Double Binomial Sum (see equation 13 in [5]). 


4. Faulhaber’s formula and finite differences 

In this section we try to go into details about our identities m , (ED , m ■ 
On the page 9 of [21] we can find an identities 

n= (i) 

n 3 = 6(”r) + (") 

n 5 = 120 ("+ 2 ) +30("J 1 ) + (?) 

For example, consider a first order finite difference applying above identities, we 
have 

= O 

in 3 = 6(“+‘) + (;) 

An 5 = 120(”+ 2 ) + 30("+') + (;) 
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These identities may have the view 


(4.1) 


E k J(n,t)n(2 n „* + TLh 


m = odd 
m = even 


where J(n, k) is defined by same identity. Particularly, coefficients J(n,k ) are 
related to what Riordan ([22], page 213) has called central factorial numbers of the 
second kind. The f-order finite difference of monomial n m is 


(4.2) 


J(n, k) ( 2n ££ t * 2fc ), m = odd 
Ek J ( n > fc)«( 2m n +i-7-2fc) > m = even 
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6. Conclusion 


In this paper particular pattern, that is triangle A287326 in OEIS, which shows 
the expansion of perfect cube n as row sum over 0 < k < n — lis generalized for 
each odd power 2 m +1, m = 0,1, 2 ..., we have reached a result 

n— 1 m 

Y Y A m j(n - k) 3 k° =n 2m+1 , m = 0,1,2,..., 
k =0 j =0 


where A m , 7 - is defined by definition (2.12) The coefficient M\(ji,k) is defined by 


definition (1.18) and generalized to M m (n,k), m > 1 at Isection 21 Properties of 


Mi(n,k) and M m (n,k) are shown in (11.201) and Isubsection 2.T1 respectively. An 
analog of MacMillan Double Binomial Sum [5] is shown in Isection 3l Relation be¬ 
tween Faulhaber’s sum ^ n m and finite differences of power are shown in Isection 41 
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